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<Summary> Shape representation plays an important role to preserve shape geometry information of an object. The 
majority of existing shape representation methods used sphere to approximate 3D object. Due to the divergent nature of the 
sphere the fitting between the sphere and the object is not tight and the mapping between sphere and object is non-compact. 
Hence the resultant spherical shape representation becomes less informative. In this paper, we propose a shape 
representation method based on an ellipsoid to represent a 3D object. The ellipsoid is convergence in nature and is better 
fitted to the body of the shape rather than spherical surface. The proposed shape representation method using ellipsoid is 
able to generate informative shape sampling function capturing detail information smoothly of an object. Our proposed 
method decreases the approximation error than the conventional spherical based method. Based on that ellipsoid shape 
representation, a Daubechies4 wavelet shape descriptor is also presented. The proposed shape descriptor provides better 
discrimination. Moreover, the approaches minimize the size of shape sampling and storage requirement. 
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1. Introduction 

Human vision can percept and define any shapes, 
images or objects especially three dimensional (3D) world 
using eyes and brain. The human learns to recognize the 
world from beginning of life acquiring knowledge of the 
world. In computer vision, it is very difficult to recognize 
objects as humans. The learning process to acquire 
knowledge may be more useful in the case of video 
surveillance system, robots and etc. Shape is one of the 
most useful and important properties to describe an object. 
Shape visualization of an object is one of the most 
important tasks and it is tolerable to humans, even parts of 
shape are not visible. Thus a challenging issue in computer 
vision is an efficient shape representation. Needlessly to 
say, we live in a 3D world. The recent development in 
camera technology such as 3D scanner and stereo vision 
leads to increase datasets of 3D model worldwide for 3D 
shape representation & recognition and gets greater 
attention for research. 3D shape representation and 
recognition can benefit many applications such as 3D 
search engine, robotics, biometric authentication, 
CAD/CAM industry, medical imaging, molecular biology 
and computer vision. 

The goal of this research is to create a mapping 
between human and machine vision by considering the 
efficiency of 3D shape representation to store their 
attribute in a highly visible manner. In this paper, we 
propose an efficient ellipsoid based shape representation 
method approximating the 3D object into an ellipsoid. The 
approximation error measured using Hausdroff distance 
using Princeton Shape Benchmark (PSB) 1} shows less 
error than the other conventional shape representation 
method based on the sphere which is used for generic 3D 
object recognition. Additionally, a robust shape descriptor 


containing Daubechies4 (Db4) wavelet coefficients 
features using the ellipsoid shape represented object is also 
proposed. 

The rest of the paper is organized as follows. In 
section 2, we present a brief review on existing shape 
representation and also state their limitations. In section 3, 
we investigate the problem of non-smooth spherical shape 
representation method used in the existing methods with 
examples of 3D shapes. In section 4, we demonstrate how 
an ellipsoid avoids the non-smooth shape sampling 
problem as a 3D shape representation method. Section 4 
also introduces an informative 3D shape representation 
method capturing detail information. Section 5 presents a 
shape descriptor which is based on Daubechies4 (Db4) 
wavelet coefficients. In section 6, we provide the 
comparison results of the proposed method against the 
existing method measuring Hausdroff based shape 
approximation error. We evaluate the descriptor in section 
7. In section 8, evaluation results are analyzed. Finally, we 
conclude the paper and state future work in section 9. 

2. Related Work 

In this section, we provide a brief review of existing 
shape representation methods which is employed for 3D 
object recognition purpose. Survey paper can find this 
topic in Ref. 2)-4). According to Loncaric 3) , shape 
representation schemes are used to preserve important 
characteristics. Shape representation schemes are useful 
since it visualize the shape accordingly. 

Conventionally, there are two categories for shape 
representation: local and global methods. Local methods 
attempt to represent objects as a set of faces or edges such 
as distance between faces or the angle between faces 
normal. Early local methods are used to handle polyhedral 
objects, but not always suitable for curved objects, 




including Ref. 5) and Ref 6). Vranic and Saupe 7) also 
presented an approach for 3D model retrieval by finding 
the number of polygons, vertices, surface area and etc. 
Earlier local methods were limited to preserving the 
object’s attribute where the features are extracted directly. 
Moreover, local methods are noise sensitive and limited to 
their polygon connectivity. 

Global methods are meant for mapping an object to a 
specified coordinate system in a predefined domain as an 
implicit function such as spherical function. Researchers 
such as Saupe&Vranic 8) , Funkhouser et al. 9) , Kazhdan et 
al. 10) , Laga et al. n) and Lmaati et al. 12) used sphere as 
shape primitive to represent a 3D shape. Afterwards, they 
employed spherical shape function to estimate shape 
descriptor based on spherical harmonics, moments, 
wavelets and spherical helices which were used for generic 
3D object recognition purpose. 

Although the shape function mapping methods on a 
sphere are popular means, they are limited to ‘fitting’ an 
object. The sphere is divergent in their nature and thus the 
representation is not closure to the body of the shape. 
Moreover, Kajya 13) mentioned that the spherical surface 
does not contain snugly. A preliminary work using an 
ellipsoid based shape representation and Haar wavelet 
based descriptor 14) is presented in our previous research. 

In the literature there are also techniques which used 
multiple numbers of primitives such as Ref. 15), 16) as a 
tree structure shape function for shape representation. In 
another method Mademilles et al. 17) proposed ellipsoid 
harmonics as a shape approximation function. In their 
method, they used 24 numbers of ellipsoids and their 
intersections with the object are regarded as a shape 
function. Then they extracted ellipsoid harmonics 
coefficients as a feature descriptor from those shape 
functions. Although those type of methods decrease 
approximation error but the complexity is also increases 
due to the requirements of number of primitives. Those 
methods may not be suitable to retrieve 3D object from a 
large number of databases. In addition, Sona 18) reported 
that numerical computation problem and unstable behavior 
of series expansion also occur in case of finding ellipsoidal 
harmonic coefficients. Moreover, extra collision handling 
method 19) is also required to avoid the collision among the 
primitives. 

3. Non-Smooth Sampling Problem on Spherical 
Shape Function 

The common and popular shape function is defined on 
the unit sphere measured on angles 6, 0 < 0 <n and cp, 0 < 
(p < 2tz. The spherical extent function (EXT) 8) / (6 , cp) is 
regarded as the shape representation function. Commonly 
Ray casting method 2 20) is used to sample the object by 
casting rays from the object’s centre in the radial direction 
(0, cp). A Ray R(t) is defined as 

R(t ) = 0 + tD (1) 

Where, O is the origin of the ray, t (>0) is the distance 
from the ray origin and D is the direction. Rays are passed 


in the radial direction and ray-triangle intersection points 
are computed. Triangle barycentric coordinates 2 ^ are used 
to determine whether the intersection point belongs to 
inside the triangle. If there is more than one intersection 
points the farthest point is taken into account and if 
intersection does not occur, then the center of mass is 
regarded as the point. Shape function is computed by 
obtaining the points at the different location. Elere, we 
analyzed the spherical shape representation function using 
the grid size 64*64 which is a popular representation. 
Figure 1(d) shows an illustration of a spherical based 
representation of examples of shark and airplane models. 
The models are collected from PSB. 

The non-uniform problem of spherical shape 
representation illustrated in Fig.l (d) may be occurring due 
to non compact fitting of the sphere to the object. Existing 
spherical mapping methods are not able to preserve the 
attribute of the shape well due to the divergent nature of 
sphere, and thus its lack of uniformity and compact shape 
representation. Spherical mapping generate a non-uniform 
representation in a divergent region of a spherical surface. 
Hence the distinguish parts of shapes for example tail and 
head part of the airplane and shark model can not provide 
a smooth shape representation which is shown in Fig.l (d). 
The resultant spherical sample shows that head, fins and 
tail parts are poorly sampled. The term ‘poorly’ means that 
the sampling has not occurred properly where the original 
shape information can not extract substantially. The non- 
uniform shape function means that central part of the 
model is well sampled whereas outmost part of the object 
provided poorly sampled. Hence this representation based 
on sphere provides less informative representation. This 
type of sampling causes loss of characteristics information 
which will not be able to provide informative 
representation of any 3D model for feature extraction and 
recognition purpose. Generally in machine vision, an 
object is more distinguishable to their outer part with the 
different objects. The outer region refers to the part or 
region of a model which is away from the centre of a 
model, for example fin, head and tail of a fish model. 
Furthermore, non-uniform shape representation may also 
lead to increase the approximation error in object 
recognition system. 

4. Ellipsoid Based Shape Representation 

To avoid the aforementioned problem explained in the 
section 3, we propose an ellipsoid based compact 
representation to generate an informative and uniform 
shape function of 3D object. The ellipsoid is convergence 
in nature and better fitted to approximate the object. The 
proposed method using ellipsoid can be used to employ 
with enriched shape information to lower the shape 
approximation error. Equation 2 shows the definition of 
an ellipsoid and the coordinate function can be used to 
compute an ellipsoid. In the Cartesian coordinate system 
an ellipsoid can be represented by 
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(c) original objects (d) spherical represented objects (e) ellipsoid represented objects 

Fig. 1 (a) spherical surface, (b) ellipsoid surface, (c) original objects, (d) spherical representation which causes non-smooth 
and non-uniform approximation, and (e) the proposed ellipsoid based objects providing smooth shape representation 


The parameters a , b and c are the semimajor and 
semiminor axis of an ellipsoid. Once the ellipsoid surface 
is computed ray-triangle intersection algorithm 21} is 
applied and the set of corresponding points on the surface 
are obtained f(6,cp), E v = [X (e>(p) Y (e>(p) Z ( e !(p) \. The following 
section 4.1 describes how to compute ellipsoid shape 
fitting. The ellipsoid is computed using the equal grid size 
64*64 providing 4096 sampled points as the spherical 
representation which has shown in Fig.l. In Fig. 1(e) we 
show ellipsoid based reconstructed objects. PCA 22) is used 
to normalize 3D models and the centre of the model is 
computed as follows: 

Consider that there is a set of vertices V of a model 
where, 

V = {Vi, v 2 , ... VnW = [x iyi zj) e S 3 (3) 

So the centre point V c (x c , y c , z c ) is defined as 

n n n 

X ‘ = 1 nY J X '' y ‘ = l nY, y <' Z ‘ = 1 n I z ‘ (4) 

i=1 i-1 i-1 

Sampling with more points may improve the quality 
of approximation. Below three reasons are discussed for 
choosing 64 as the used sampled dimension: 


(i) Firstly, it is observed that, if the shape is sampled 
using 32 then the sampling provides a poor 
approximation for an example of an ant model (from 
PSB dataset) using the ellipsoid approximated shape 
function. It is also noticed that shape sampling with 
64 and 128 did not show significant difference in the 
ant model which will make us to consider increasing 
the sampling rate. Increasing the sampling rate may 
provide a better solution for more complex objects. 

(ii) Secondly, space complexity is higher to the 
exponential, if the sampling rate will be regarded as 
128. 

(iii) Thirdly, time complexity also increases when the size 
of the vertices of models increases. For example a car 
model containing 308 numbers of vertices took 5.9 
seconds to sample the object when the sampling rates 
were 64. The computation time took 22.2 seconds, if 
the sample size were selected as 128 for the same 
model. Figure 2 shows the extraction time verses 
sampling with various numbers of vertices of 
different models. 

4.1 Ellipsoid surface computation 

The parametric equation of an ellipsoid can be 
defined 23) as 

k = a cos u sin v 
l = b sinu sin v 
m = c cos v 


( 5 ) 
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Fig. 2 Extraction time increases, according to sampling numbers 



(a) sphere based 


(b) ellipsoid based 


Fig.3 The difference between proposed ellipsoidal and spherical representation, points are distributed after 
intersections to the primitive 

















where u E [0, 2i r], v E [0, 7r] and a , b and c and V a, b & c 
^ 0, are semi major axes. 

If the lengths of two axes of an ellipsoid are same, it is 
spheroid, while a > c called oblate spheroid and when c > 
a it is prolate spheroid 23) . A sphere is a special case of 
ellipsoid when a = b = c. The ellipsoid surface E (X e , Y e , 
Z e ) is computed using the equation (6) and (7). 

X e = a*x, Y e = b*y , Z e = c*z (6) 

where, 

x= r sin 0 cos cp 

y = r sin 6 sin cp (7) 

z = r cos 6 

with a radius r and (0<cp< 2n) and ( 0<0<n ) 

The ellipsoid axis parameter such as a , b and c are set 
into 1.5, 1 and 0.5 respectively, since it is observed that 
this combination gave the best results. From the 
observation, the proposed mapping algorithm must set the 
axis parameter as a > b > c, otherwise poor sampling 
occur. 

4.2 Ellipsoid as a uniform shape function 

Ellipsoid shape representation avoids the non-smooth 
sampling problem occurred in the case of spherical 
sampling, which is illustrated in Fig.l and Fig.3. The 
existing spherical surface sampling produces irregular 
sampling. The term ‘irregular’ means denser sampling in 
the center part and poor sampling in the upper and lower 
comer parts of a 3D model. It is known that sphere is 
divergence nature due to pole exists and ellipsoid is 
convergence in nature due to semi major and semi minor 
axis presents in the ellipsoid. Moreover ellipsoid does not 
consist any pole. Figure 3 shows a difference between 
spherical and ellipsoidal surface (shown in points only) 
mapping of two different model ‘shark’ and ‘airplane’. 
The points represent the distribution after spherical and 
ellipsoidal sampling of those models. As the sphere does 
not fit well to the object and hence the distance between 
the intersection point and ray origin becomes higher. 
Hence the points near to the poles of the model are going 
far from points to the origin of the surface. Hence the 
chances of ray object intersection decrease according to 
the criteria of shape approximation 24) and those points (far 
points) are poorly sampled. The cases while point’s 
intersection does not occur, center points are taken as an 
intersection according to ray tracing procedure 20) . Hence 
the center point of the spherical mapping goes denser due 
to the above occurrences. Due to the convergent nature of 
the ellipsoid, the points away from the center of the model 
are coming nearer to the origin of the ellipsoid surface. 
Moreover, the detail parts of the 3D model come nearer to 
the center since the distance between the intersection point 
and the ray origin becomes nearer due to close fitting. 
Eventually, the detail parts sampled uniformly similar to 
the center part of a model using the proposed algorithm 
than sphere based sampling. This generates the 
informative, smooth, regular and compact shape function 
using the proposed ellipsoid mapping. From the 
observation, the intersection points are denser in the 



(a) (b) (c) 


Fig.4 Render images of (a) an ant and (b) a flying 

bird and (c) a tree model applying the proposed 
Ellipsoid Approximation method 



ellipsoid represented image decomposed image after 
of‘fighterjef model, image level 1, image size 32*32 
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decomposed image after 
level 2 ( 16 * 16 ) 


Fig.5 A 2 -level Db4 wavelet transform decomposition 
of an ellipsoid mapped image of a ‘ fighterjef 
model 

distinguish parts of 3D models in the case of the ellipsoid 
representation rather than spherical representation. Hence 
the distinguish parts (outmost) which are away from centre 
part of any object can be preserved better with the 
proposed representation. 

The proposed method not only approximates shape 
better but also minimizes the size of the shape sampling. 
In addition, the dimension of the descriptor will also 
minimize if the sampling dimension minimize. In general, 
ray casting method is slow and the time complexity also 
increases using ellipsoid since the sampling requires extra 
operation due to the semimajor and semiminor axis of 
ellipsoid. Still the performance is not degraded. 

5. Wavelet Descriptor Extraction 

A feature descriptor is made by applying Daubechies4 
(Db4) 25) wavelet transform to describe a 3D shape on the 
ellipsoid represented object. Db4 wavelet coefficients are 
extracted from level 2. The feature vector contains 512 
numbers of detail and approximation coefficients. 







Before extracting descriptor, we have unfolded the 
ellipsoid represented shape function into 2D image since 
the resultant ellipsoid shape function is high dimensional 
and prone to noise effective. The obtained image is used to 
employ feature descriptor. Examples of rendered images 
are shown in the Fig.4. The Fig.5 shows the decomposed 
image of a 2-level Db4 wavelet transform of an ellipsoid 
represented object model. We used Euclidean distance for 
similarity matching of two feature descriptors. 

6. Approximation Error Measurement 

This section presents the results on shape 
approximation error measurement and comparison result 
between existing spherical and proposed ellipsoid based 
shape representation. We used Hausdroff distance 26) to 
measure the approximation error. Hausdroff distance has 
been used for measuring surfaces error in many 
applications such as mesh simplification, geometry 
compression and surface distortion 19) . Hausdroff distance 
d s betweenthe original surface X and the approximated 
surface Y can be defined as: 

d s (X, Y)= max [d(X, Y), d(Y,X )] (8) 

Where the relation d (X, Y) =™exd (p, Y) is between a 
point p belonging to a surface X and a surface Y is defined 


as d(p, Y)= II (p — p') II2- 

We compare the approximation error of the proposed 
ellipsoid based shape representation against spherical 
representation with the example of various fighter Jet, 
hand gun, fish and bird models of Princeton Shape 
Benchmark (PSB) models. Those models are chosen 
randomly from the database. A class contains various 
types of models for example, fighter Jet class contains 
different types of fighter Jet models with different 
orientation etc. In Fig.6 (a), the horizontal axis represents 
various fighter Jet modelsfrom fighter Jet class which is 
taken randomly from Princeton Shape Benchmark, and so 
on for the Fig. 6 (b)-(d). 

From the observation, the shape approximation error 
decreases using the proposed ellipsoid representation 
method. For example, Fig.6 (a) shows that the error is as 
low as 30% of a fighter model using ellipsoidal rather than 
spherical representation. It is to believe that due to the 
compact representation of ellipsoid tothe shape, the 
approximation error is smaller than the existing sphere 
based approximation. 

7. Descriptor Evaluation 

In this section, we have shown the performance and 
comparison results of the proposed shape descriptor based 
on the ellipsoid shape representation method. The 
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Fig. 6 Comparison of approximation errors of the proposed ellipsoid shape representation scheme 

against spherical representation of various examples: (a) fighter Jet, (b) handgun, (c) fish and 
(d) bird models 





















































computational time 0.070 seconds require for building a 
descriptor of the ellipsoid represented image. We used 
standard evaluation tools such as Nearest Neighbor (NN), 
First Tier (FT) and Second Tier (ST) to ensure 
effectiveness of the proposed descriptor. Table 1 
summarized the results on the PSB base test classification 
which contains 907 models of 92 classes in the database. 
The test in alternate classification, coarser classification is 
also examined which contains 6 classes and is presented in 
the Table 2. The effectiveness of the proposed descriptor, 
ellipsoid wavelet descriptor (EWD) is measured 
comparing with other well known spherical based methods 
such as spherical wavelet descriptor (SWD) U) ’ spherical 
harmonic descriptor (SHD) 10) , spherical extent function 
(EXT) 8) , discrete cosine transform (DCT) 12) . The Table 
3 lists the comparison performance on PSB base test 
classification. The comparison results between spherical 
and proposed ellipsoid based on coarser classification are 
shown in Fig.7 (a), (b) and (c) using NN, FT and ST 
respectively. 

8. Analysis for Results 

This section analyzes the experimental results of the 
proposed descriptor EWD based on discrimination, storage 
requirement and time requirement. 

8.1 Discrimination analysis 

The overall performances of the proposed EWD have 
shown an improvement results from all other descriptors 
using all the measurements. It can believe that the overall 
performance of the proposed method increases since the 
shape representation method using ellipsoid provides 
smooth and informative representation which increases the 
discrimination attribute and outperforms in shape 
matching. Hence the performance increases for the 
majority of the classes. Additionally, the Db4 wavelet 
features also provide effectiveness of the performances 
and robustness of the descriptor, due to its compact 
support. More specifically, the result shows that the 
proposed descriptor attains good performance in the 
vehicle, buildings, plant and few of animal classes other 
than household class due to their more compact shape 
attributes. 

8.2 Analysis based on storage requirement 

The Table 4 shows the comparison of storage 
requirements for shape sampling function and descriptors. 
This comparison shows that the proposed shape 
representation method used fewer shape sampling size 
than the spherical based approach although the proposed 
method is able to provide good informative representation. 
Moreover, it shows that proposed descriptor EWD has 
fewer dimensions against all the descriptors except EXT. 
The proposed descriptor EWD outperforms with notable 
performance than those which have presented above 
although the descriptor is able to minimize the storage size. 
The proposed descriptor in this research is able to take the 
benefit of being lessexpensive for real time applications. 


Table 1 Performance of the proposed descriptor EWD 
of PSB datasets base test classification 


Descriptor 

NN 

FT 

ST 

Proposed 

(EWD) 

61.1% 

33.1% 

41.5% 


Table 2 Evaluation results of the proposed ellipsoid descriptor 
in PSB coarser test classification (6 classes) 


Query Classes 

Nearest 

neighbor 

First tier 

Second tier 

Vehicles 

78.36% 

44.89% 

69.38% 

Animals 

67.74% 

32.25% 

49.03% 

Buildings 

59.57% 

12.77% 

19.15% 

Furniture 

51.06% 

11.70% 

23.4% 

Plants 

58.33% 

10% 

16.67% 

Household 

34.59% 

18.37% 

25.41% 


Table 3 Performance comparison of the proposed EWD using 
NN, FT & ST on PSB test base classification 


Descriptor 

Nearest 

neighbor 

First tier 

Second tier 

Proposed EWD 

61.10% 

34.4% 

41.0% 

SWD 

46.90% 

31.4% 

39.7% 

SHD 

55.60% 

30.90% 

41.1% 

EXT 

54.90% 

28.60% 

37.9% 

DCT 

53.50% 

26.5% 

36.4% 


Table 4 Comparison of the ellipsoid & spherical approximation 
methods and descriptors according to size 


Method 

Surface grids 
size (S) 

Sampling points 
(S*S) 

Descriptor 

storage 

requirement (in 
byte) 

Proposed EWD 

64 

4096 

512 

SWC 

128 

16384 

512 

EXT 

128 

16384 

136 

DCT 

190 

4096 1 

1368 


Noted that the sampling points with spherical helix relations are 
based on helix parameter q, e.g. if q =190, points are 4096, more 
illustrated in the paper 12) 


8.3 Time measurement analysis 

Figure 7 (d) shows the time taken for shape 
representation method of the proposed ellipsoid 
representation and also compare against the spherical 
representation. It shows that the proposed ellipsoid based 
approach takes more time than spherical sampling which 
can be regarded as a limitation. The reason is increasing 
the extra operation for the semi axis exist. Still the 
performance is not degraded. 
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Comparison evaluation of EWD against SWC using ST in 
coarser2 classes 
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Fig. 7 Performance comparison of EWD against SWC using (a) NN measure (b) FT measure (c) ST measure on PSB test coarser 
classification, and (d) comparisons of time taken for sampling using sphere and ellipsoid based mapping method, where 
surface grid size is 64*64 for both of the mapping 


9. Conclusion 

In this paper, we propose an efficient uniform and 
compact ellipsoid based 3D shape representation method 
and based on that method, a shape descriptor of Db4 
wavelet coefficients is also presented. The shape 
representation method provides an informative 
representation which is very useful and efficient for 3D 
shape descriptor extraction purpose in general kinds of 
shape. The non-smooth sampling in the case of spherical 
mapping is avoided by introducing the idea of proposed 
ellipsoid shape representation. Moreover, the ellipsoid 
shape represented method lowers the approximation error. 
Additionally, the Db4 based wavelet descriptor provides a 
compact and robust descriptor for shape description. 
Specifically, the proposed descriptor performs well in case 
of vehicles, buildings, plant and few of animal classes of 
the PSB database. The advantage of the proposed shape 
representation as well as the descriptor has fewer 
dimensions in size. The limitation of the proposed 
ellipsoid mapping scheme is that sampling takes more time 


than spherical since ellipsoid sampling needs some extra 
operation for ellipsoid coordinate axis. 

Moreover, we are considering this limitation as a 
future work to apply a grid accelerating structure method 
on the ray tracing approach currently used in this paper. A 
more compact and fast like uniform grid structure or 
hashing grid structure or any other suitable tree structure 
might be applied to optimize the effectiveness of this 
proposed method. Using these accelerated grid structured 
methods the testing of the intersections between the object 
and the ellipsoid will be reduced to those only, where the 
intersection occur. 
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